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SOME STUDIES ON THE NONLINEAR DYNAMIC RESPONSE 

j 

OF SHELL-TYPE STRUCTURES 
By David A. Evensen and Robert E. Fulton 
NASA Langley Research Center 

INTRODUCTION 

Recent aerospace applications involving thin-walled structures have neces- 
sitated studies of the nonlinear vibration and response of thin shells. The 
results of some current research on nonlinear dynamics of thin shells are sum- 
marized in the present paper. Two related problem areas are outlined; both 
involve geometric nonlinearities and thin shells of revolution. 

The first problem considered is the forced, nonlinear flexural vibrations 
of thin-walled circular cylinders. Approximate solutions are obtained which 
show that vibrations involving a single bending mode or two coupled bending 
modes can occur. The analysis exhibits several features that are character- 
istic of nonlinear flexural vibrations of axisymmetric elastic structures in 
general . 

The second problem deals with the dynamic axisymmetric snap-through buckling 
of shallow conical and spherical shells subjected to uniformly distributed impul- 
sive pressures. Approximate buckling pressures are obtained for a variety of 
boundary conditions for both types of shells and comparisons are made between 
conical and spherical caps having similar geometries. 

SYMBOLS 

An('t),B n (t) generalized coordinates 
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average values (over one period) of the vibration amplitudes 
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(see eqs. (7) and (9)) 
external radius of conical shell (fig. 7 a ) 


bending stiffness. 


Eh^ 

12(1 - b 2 ) 


Young 1 s modulus 
stress function 

nondimensional amplitude of the applied loading 
shell thickness 

center rise of the shell (figs. 7 a and 7 b) 


initial impulse parameter, 


^bSd - b 2 ) 

Eph^H 2 


initial impulse per unit area 

radial and circumferential changes in curvatures 
length of the cylinder 

number of axial half-waves and circumferential waves, respectively 
radial loading applied to the surface of the cylinder 
horizontal coordinate of conical shell (fig. 7 a ) 
radius of the shell 
time 

displacements of a point on the median surface of the shell (see 
figs. 1 , 7 a, and 7 b) 

displacement constant, see equation ( ll) 
shell coordinates (see fig. l) 
nondimensional strain energy of the shell 
colatitude coordinate for spherical shell (fig. Tb) 



half opening angle for spherical shell (fig. 7b) 




nonlinearity parameters 


« 1^2 


radial and circumferential extensional strains 


£c( t )>£s( t ) nondimensional generalized coordinates associated with and 

Bjj, respectively 


shell parameter, \ S ^8( X - 


Poisson’s ratio 


aspect ratio of the vibration mode, 


mass density 


radial and circumferential stresses 


nondimensional time , cngt 


vibration frequency 


calculated linear vibration frequency defined, by 
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e(t 2 + l) 2 lV /2 


PR 2 (|2 + i) 2 12(1 - n2) 


experimental linear vibration frequency 


nondimensional frequency, — 


biharmonic operator, [ + 2 + 

\Sx^ chc 2 c)y 2 Sy^y 
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NONLINEAR FLEXURAL VIBRATIONS OF THIN -WALLED CIRCULAR CYLINDERS 

The nonlinear flexural vibrations of thin-walled, circular cylinders are 
analyzed by assuming two vibration modes and applying Galerkin's procedure on 
the equations of motion. The assumed shape for the radial deflection w and 
the related stress function F are chosen such that they approximately satisfy 
the boundary conditions of a cylinder having "freely-supported" ends. In-plane 
inertia terms are neglected, and a linear stress-strain law is assumed. The 
nonlinearities in the problem arise from including nonlinear rotation terms in 
the strain-displacement relations. 


Governing Equations 

Using the well-known approximations of Donnell's shallow-shell theory as 
exemplified in reference 3, the equations of motion of a thin-walled circular 
cylinder can be combined to give 


dAw + ph 

St 2 


q + ±£l + 

R 5x 2 


<£f cfw 

by 2 dx 2 


2 d 2 w + d 2 w 

dx by dx by £ x 2 $y 2 


( 1 ) 


and 



i a 2 w + 
R Sx 2 



d 2 w d 2 w 
bx 2 by 2 


( 2 ) 


where w is the radial deflection and F is the usual stress function. (The 
coordinate system and shell geometry are shown in fig. 1.) 

Approximate solutions to equations (l) and (2) were obtained by using the 
following two-mode approximation for the radial deflection: 
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w(x,y,t) = jl n (t)cos + Bn(t)sin S^sin ^ 

+ ^lA n 2 (t) + Bn 2 (t)ls: 




Here cos — sin 2J2E. and sin — sin 222 are the linear vibration modes of 
R L R L 

the cylinder, and equation ( 3 ) is limited to n > 2 since only flexural motions 
are being considered. The bracketed term in ( 3 ) is included to satisfy the 
periodic continuity condition on the circumferential displacement v. 

Substitution of equation ( 3 ) into the compatibility equation (2) allowed 
the latter to be solved for F. The expressions for w and F were then 
examined and found to satisfy the following boundary conditions: 

(a) The displacements u, v, and w, and their derivatives satisfy perio- 
dicity conditions of the form 


v(x,y,t) = v(x,y + 2*R,t) 

(b) The radial displacement w goes to zero at the ends of the cylinder, 
i . e . , at x = 0 and x = L. 

(c) To a first approximation, the in-plane stress-resultant N x , the 
moment -re suit ant M x , and the tangential displacement v all vanish at x = 0 
and x = L. In other words, the coefficients of the linear terms in the expres- 
sions for N x , Mx, and v go to zero at the ends of the cylinder, but the 

O O 

nonlinear terms involving An", A n Bn, and Bn do not vanish there. These 
end conditions are similar to the boundary conditions for a cylinder which has 
freely-supported ends. 

The applied loading, q(x,y,t), was chosen such that only one mode is 
directly excited: 


5 


/ , x rxY m jtx , 

q(x,y,t) = Qmn cos y sin — cos “* 


W 


Finally, the expressions for w, F, and q were substituted into equation (l) 
and a Galerkin procedure was used to obtain two nonlinear differential equations 
for the modal amplitudes A n and Bn. In nondimen si onal form, these coupled 
equations are 




dT £ 


+ «« + f 


d 2 £ /d 5 \ d d £ /dt 

tc + ( — ) + + — 

° dT 2 \dx / Ss dT 2 \dr 


^c(t>c 2 + ^s 2 ) + e2 ^ c (^c 2 + ^s 2 ) 2 = G mn ^ Or 


(5a) 


and 


dr^ 


r + t 

T o Sc 


+ ls-J + ?c 
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(5b) 


where the nondimensional variables are 
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- h 
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^s h 


T = 0 ) g t 
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mu q q 
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the nondimensional frequency is 


ft = 


60 

OOg 


and the nonlinearity parameters are 


e = 


[nfh 

V R ) 
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and 


5 = 


5 | 

16 r 


[(I 2 + l) 2 (9I 2 + l) 2 ] 


e(l 2 + 1)‘ 


(|2 + l ) 2 12(1 - u 2 ) 


In these expressions, the parameter | represents the aspect ratio of the 
vibration mode: 


^ _ itR/n _ Circumferential wavelength 
L/m Axial wavelength 


(6) 


Examination of the expressions for 7, 8, e, and i in equations (5) 

and (6) leads to the following observations: 

(a) As the length of the cylinder tends to infinity, the parameters £, 

7, and 5 all tend to zero, and equations (5) approach the previous results 
for rings (ref. l) . 

(b) Each of the nonlinear terms in equation (5) is multiplied by e; con- 
sequently, e can be viewed as the basic nonlinearity parameter in the problem. 
Linear vibrations occur for e = 0, and increasing € makes the vibrations 
increasingly nonlinear. 


Approximate Solutions by the Method of Averaging 
Equations (5) can be solved approximately by use of the method of Krylov- 
Bogoliubov, often called "The Method of Averaging." Such solutions will be 
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presented for vibrations involving (a) only the driven mode and (b) both the 
driven mode and its companion mode. In the discussion which follows, the term 
"driven mode" relates to £ c (t) and cos ny/R sin nmx/L; the term "companion 
mode" refers to £ S ( T ) and sin ny/R sin mitx/L. 

Response of a single bending mode .- Since the applied loading (eq. (4)) 
directly drives only one mode of the cylinder, a possible solution to equa- 
tions ( 5 ) involves the response of only the driven mode. Application of the 
method of averaging for this case gives 


£ c (t) = A cos £2t[ 
C s (t) = 0 


where A can be computed from 


(7) 


(1 - fi 2 )A - n2 a 3 _ 5e7 a 3 + 5 e 2 5A 5 = G ( 8 ) 

16 4 8 

Equation ( 8 ) can be used to compute the variation of A with 0, for given 
values of e, 7 , 6 , and G^. When is nonzero, the forced vibration 

response of a single mode is obtained; the case of free nonlinear vibrations 
results when G mn is put equal to zero. 

Stability of the one -mode response .- The stability of the preceding solu- 
tion was investigated by perturbing £ (t) and £ s (r). A study of the 
resulting Mathieu-Hill equations indicated that within terms of order e 2 
(l) Perturbations of £ c are unstable within the area bounded by 

1 • 9 j fk + i) K a K 1 • ( r * i) 
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(2) Perturbations of £ s are unstable within the region 


1 - 


JeA 2 

8 


( r + j) <a 


<-#(H 


a = 


(5) Both types of perturbations are unstable in narrow regions near 

1 / 2 , 1 / 5 , . . . 


The first instability region coincides with the locus of vertical tangents 
to the response curves and indicates the well-known jump phenomena. The nar- 
row areas near ft = l/2, l/5, . • ., denote possible ultraharmonic responses. 
The remaining region, (2), indicates the area in which the companion mode is 
parametrically unstable due to nonlinear coupling with the driven mode. To 
obtain adequate solutions in region (2), it is necessary to consider motions 
where both modes vibrate. 

Response of the coupled bending modes .- When £ s (t) and £ c (t) both 
oscillate, the method of averaging gives the approximate solution 


£ c (t) = A cos ftT | 
£ s (t) = B sin ftT | 

where A and B satisfy the following equations: 


(9) 


(l - ft2)A + A(l2 - I 2 ) - ■^(5A 2 + B 2 ) + ^ a(5A^ + 2A 2 ! 2 + B^) = 

16 4 o \ 

(10a) 

(i - n2)s + b(a^ - #) - iZ2(s s 2 ♦ a 2 ) + £5. 1(5/ + as 2 ! 2 + a 4 ) = o 

(10b) 

If B = 0, equations ( 9 ) and (10) revert to the single mode case discussed 
previously. When B is not zero, the variation of A and B with ft can be 
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computed by solving equations (10a) and (lOb) simultaneously, for given values 
of e, 7, 6, and G^. 

DISCUSSION OF THE CYLINDER RESULTS 


Figure 2 illustrates the response curves for both free and forced vibra- 
tions of a single bending mode. The dashed curve represents free nonlinear 
vibrations and was obtained from equation (8) with Omn put equal to zero. 

The forced vibration response is given by the solid lines, which were computed 
from equation (8) for G rrm = 0.1. Both the free and forced response curves 
demonstrate a slight nonlinearity of the softening type. The value of e used 
in these calculations was 0.01, and the values of y and 5 which were used 
correspond to i = 0.1, e = 0.01, and p = 0.3. These values of | and e 
are representative of a cylinder that is vibrating in the m = 1, n = 10 mode 
and which has a length/radius of n and a radius/thickness ratio of 1000. The 
slight nonlinearity shown in figure 2 is typical of cylinders that are rela- 
tively long and thin-walled. For other modes and geometries, however, the 
nonlinearity can be much more pronounced. 

This result is indicated in figures 3 and 4, which show how variations in 
the parameters e and £ affect the nonlinearity of the vibrations. Fig- 
ure 3 shows several free vibration response curves computed from equation (8) 
for Gum = 0 and for five values of e ranging from 0 to 1.0. The solid 
lines were calculated for values of y and 5 corresponding to an aspect 
ratio | = l/2, and the dashed curves are for | = 2. Both sets of curves in 
figure 3 demonstrate that the strength of the nonlinearity is determined pri- 


marily by the parameter 



The nonlinearity is generally small for 
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vibrations involving very thin cylinders and/or low circumferential mode num- 
bers, n. Conversely, strong nonlinearities occur for the case of thick cylin- 
ders and/or high circumferential mode numbers. 

The character of the nonlinearity (i.e., whether it is softening or 
hardening) depends primarily on £ , which is the aspect ratio of the vibration 
mode. This result is illustrated in figure 4, which shows several free vibra- 
tion response curves computed from equation (8) for values of % ranging from 

0. 10 to 4.0. The solid lines show the results of the present analysis; they 
were calculated for constant values of the nonlinearity parameter e = 1.0 and 
Poisson’s ratio p = 0.5- The present calculations show that the vibrations 
are generally of the softening type when the aspect ratio | is less than 
unity. For larger values of the aspect ratio, the solid curves indicate a 
nonlinearity of the hardening type. 

For comparison purposes, the dashed curves in figure 4 illustrate Chu's 
results (ref. 4), which correspond to e = 1.0 and p = 0.318- Chu points 
out that his results possess a symmetric dependence on the aspect ratio param- 
eter, whereby his curves for | = l/2, l/4, l/8 . . . coincide with those for 
1=2, 4, 8. . ., respectively. Such a symmetric dependence on the aspect 
ratio seems to conflict with the basic geometric nonsymmetry of the problem - 

1. e., the surface of the cylinder is curved in the circumferential direction 
but not in the axial direction. 

With regard to this point, it will be noted that the present results do 
not exhibit a symmetric dependence on the aspect ratio and do not agree with 
the calcinations of ref. 4. The major reason that the previous studies 
(refs. 3 and 4) do not agree with the present work appears to be that the former 
do not satisfy the necessary geometric continuity constraint on the 
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circumferential displacement, v. (The available experiments tend to favor 
the present analysis; they will be discussed shortly.) 

The previous studies did not consider coupled mode responses involving 
both a driven mode and its companion. The solid lines shown in figure 5 
illustrate such a coupled mode response for a typical cylinder. The response 
curves were calculated from equations (10a) and (lOb) with e = 0.01, 

Gmn = 0.1, and y and 5 corresponding to | = 0.1. It is of interest to 
note that the solid response curves are analogous to the results obtained for 
nonlinear vibrations of rings and for nonlinear vibration absorbers (refs. 1 
and 5 * ) 

Along the a-b portion of the solid curves (fig. 5)> B is zero and only 
the driven mode responds. Along the segments b-c-d and e-f-g, both the driven 
mode and its companion respond with comparable amplitudes; in this case, the 
response curves are adjacent to the dashed lines, o-n. The response curves 
exhibit vertical tangents at points c, f, g, and h, and the segments c-d and 
f-g are suspected of being unstable . Along g-h-i-j, the response of the com- 
panion mode is much greater than that of the driven mode. This situation is 
analogous to the vibration absorber response in which the driven mass experi- 
ences very little motion while the absorber mass vibrates with large displace- 
ments. For ring vibrations, the segment corresponding to g-h-i-j was found to 
be unstable, and a similar result appears true for cylinders. To the left of 
point g, the vibrations apparently revert to the one mode case, with B = 0 
and A given by equation ( 8 ). For seme cylinders, this may result in a "gap" 
region, in which both solutions ( 7 ) and ( 9 ) are unstable. Such unstable gaps 
in the response have been observed in related problems, and analog computer 
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studies indicate that nonsteady vibrations with rapidly changing amplitudes 
occur in these regions. (See refs. 1 , 6 , and 7 in this regard.) 

Experimental results for the nonlinear flexural vibrations of thin-walled 
circular cylinders are rather scarce. Among the best available data is that of 
Olson (ref. 2) who obtained nonlinear forced response curves for the vibration 
of a thin cylindrical shell. Olson's experimental data is shown by the circles 
and dashed lines in figure 6 . The solid line in figure 6 is the free vibration 
response curve calculated from equation ( 8 ) for Gmn = 0 and with values of 
e , 7 , and 6 that correspond to Olson's experiment. Both theory and experi- 

ment indicate a nonlinearity of the softening type for this case, whereas the 
previous theoretical results (refs. 3 and. 4) indicate a hardening type of non- 
linearity for all cylinders. The agreement shown in figure 6 is reasonably 
good, despite the fact that the theory is for freely- supported ends whereas the 
test cylinders had end clamping. It is of interest to note that Olson also 
detected nonsteady vibrations in his experiments; these vibrations may have 
been due to coupled mode responses. 

Nonlinear' Vibrations of Related Axisymmetric Systems 

The results for thin cylindrical shells are qualitatively similar to those 
for rings; in both cases, the companion mode can be parametrically excited and 
participate in the motion. In some cases, the driven mode and its companion 
combine to produce a circumferentially travelling wave; responses of this type 
have been detected for rings (ref. l), cylindrical shells (ref. 8 ), conical 
shells (ref. 9)> and. circular plates (ref. 10). It seems reasonable to expect 
that similar motions will occur for the nonlinear vibration of other flexible 
axisymmetric bodies. This suggests the possibility of formulating a general 
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theory for the nonlinear vibration of axisymmetric systems, and the present 
work is being extended along these lines. 

DYNAMIC SNAP -THROUGH BUCKLING OF ELASTIC SHALLOW CONICAL 

AND SPHERICAL CAPS 

A second class of nonlinear problems which has been studied involves the 
axi symmetric snap-through buckling of shallow conical and spherical caps sub- 
jected to uniformly distributed impulsive loads. For this problem, it is 
assumed that the duration of the impulse is short with respect to the period 
of the fundamental flexural mode, but long with respect to the first exten- 
sional mode of the shell. The shell is assumed to behave as a single-degree- 
of -freedom system during its first cycle of oscillation and inplane inertia is 
neglected. 

The method of attack is to investigate the behavior of the shell at the 
time when its amplitude reaches a maximum. Treating the shell as an undamped 
single-degree-of -freedom system, the internal strain energy at maximum dis- 
placement is equal to the initial kinetic energy of the shell. The initial 
kinetic energy can in turn be related to the applied impulse by considering 
conservation of momentum. 

Thus, the initial impulse can be obtained as a function of the maximum 
displacement. The maximum displacement w will increase with increase in 
initial impulse I until the impulse is sufficient to cause snap -through 
buckling. At the snap-through impulse there will be a sudden increase in maxi- 
mum displacement with respect to a change in initial impulse. The condition 
for snap-through is therefore given by 
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Conical Shell 


The problem to be considered first is that of a shallow conical shell 
which deforms axi symmetrically. The results will be summarized briefly and 
further details may be obtained in reference 11. The cone has a height H, an 
external radius b, and thickness h. The radius to a point on the middle sur- 
face is r (see fig. 7 a )* Tbe shell is considered shallow in the sense that 
the rise angle of the cone and the sine of the rise angle can be replaced by 
H/b. 

The strain energy of the shell can be expressed as 


V = [ fe-. 2 + Gp 2 + 2pe-, Gp^rdr + f (ki 2 + K2 2 + ^ K l K 2) rdr 

2(1 - 4 2 )h 2 J 0 V X ' J 0 V ' 

( 11 ) 

where the factor 2rtH^Eh has been used to render the strain energy nondimen- 
sional. Here e^, eg an l K^, Kg are, respectively, the extensional strains 
and curvatures in the radial and circumferential directions. Young's modulus 
is denoted by E and Poisson's ratio by n; A is the shell parameter, which 
is defined by 


x 2 - !tai - 1*2) 


( 12 ) 


The nonlinear strains and curvatures given in equation (ll) are related to the 
displacements by 


_ du , H dw 
" d? b d? 




K-, = 


d^w 

dr2 




yr _ 1 dW 

Ks r dr 


( 13 ) 
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where u and w are the maximum deflections in the horizontal (r) and vertical 
(z) directions, respectively. 


Solutions for Various Boundary Conditions 
Consider first the solution for the conical cap which has boundary condi- 
tions at r = b such that the displacement is restrained in the radial direc- 
tion (i.e., u = 0). The shell can be either clamped or simply- supported with 
respect to the w displacement at the boundary. 

The maximum deflection can be approximated by 



(14) 


where w G and C are constants. The boundary conditions at r = b can be 
satisfied by a proper choice of C, that is. 


Clamped: 



C = 1 


Simply supported: 


Mq = -D(Kp + pK 2 ) =0, C = 1 + ~ ^ 

5 + p 


The shape given by equation (l4) corresponds to the deflection of a circular 
plate subjected to a uniform lateral load and supported at the boundary. 

Following the procedure just outlined leads to an equation such as the 
following one which holds for the clamped restrained boundary conditions 


i -i dw 

i.e., — 

1 dr 


= uj = 0 |. 
r=b r=b 


= >4 


0 . 0288 w o ^ - 0.0908 w o 5 + w Q 2 (0.0719 + 


5^5353 


(15) 
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where 

T .. rV*3(i - m£) 

Eh^H 2 


(16) 


I is the initial uniformly distributed impulse per unit area, p is the mass 
density of the shell, and value of p. = 0.3 has been used. 

The condition for snap-through is dl/dw c = 0 which gives the critical 
deflection for the clamped case as 


w Q = 1.1823 - \|o. 1^6 - 92-5926 


( 17 ) 


Note that at certain values of A the amplitude w Q becomes imaginary; 
at these points there is no real w Q corresponding to the snap-through condi- 
tion. The physical interpretation of imaginary values of w Q is that the 
shell does not exhibit a snapping phenomenon and instead the motion is smooth 
and oscillatory. The mini mum values for A for the clamped case corresponding 
to real w 0 is A^n = 5*0. Thus, for example, in the case of a cone with a 
clamped restrained boundary, snap-through buckling does not occur unless A 
exceeds 5*0 which happens when the rise of the cone H exceeds 3*78h* 

A plot of the impulse required for snap-through as a function of the shell 
parameter A is presented in figure 8. These curves are obtained from equa- 
tion (15) taking account of equation ( 17) • Similar calculations were carried 
out for all combinations of simply-supported or clamped and restrained (i.e., 
radial displacement u = 0) or unrestrained (i.e., radial stress o± = 0) bound- 
ary conditions. These are also given in figure 8. 

A study of the magnitude of critical deflection occurring at snap-through 
indicates that it is of the order of the rise of the cone. For example, when 
A = 6.0, — = 0.902 for the clamped restrained case. 
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Spherical Cap Results 

Results were also obtained for the snap-through buckling of a spherical 
cap having the same boundary conditions and loading as the cone. The appropri- 
ate strain-displacement relations were taken as 





k 2 = I^ 

^ ^2 a da I 

Here a is the opening angle variable and R is the radius of the sphere 
(fig. 7b) . The quantities u and w are the tangential and normal displace- 
ments of the shell. 

The displacement w is taken in the same form as equation (l4) and the 
procedure for the spherical cap is directly analogous to that for the conical 
cap. A plot of the impulse required for snap-through as a function of shell 
parameter A is given in figure 9 for boundary conditions comparable to those 
of the cone, i.e., simply-supported and clamped, restrained and unrestrained. 
The results for the clamped restrained spherical cap agree identically with 
those given in reference 12. 


Discussion of Results 

It is interesting to note in figure 8 that for relatively large values of 
A the strength of the simply- supported restrained conical shell is greater 
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than the comparable fixed case. On the other hand for small values of ~h where 
the rise of the cone is quite small, the reverse is true. 

The first effect has been noted previously by Humphreys and Bodner 
(ref. 12) for dynamic buckling of a long cylindrical panel. They found that 
for all ranges of the curvature parameter the simply-supported panel exhibited 
a greater resistance to snap-through than did clamped panels. 

The results for the spherical cap in figure 9 follow the same trend with 
the simply-supported shell always being stronger with respect to dynamic 
buckling than the clamped shell. These results emphasize the importance of 
boundary conditions in assessing the dynamic strength of this class of shells. 

It is also of interest to compare the behavior of a shallow conical shell 
with that of a shallow spherical shell of comparable geometiy and boundary con- 
ditions. Since a shallow sphere cannot be compared directly with a shallow 
cone, two comparisons were made, first with equal rises and the second with 
equal edge slopes. Both results are plotted in figure 10 for the restrained 
boundary conditions where the spherical cap parameters have been used as a 
basis for the plots. In general, figure 10 shows that for a conical or spheri- 
cal shell of the same rise H, the stronger structure to resist dynamic buckling 
is the sphere. Figure 10 also indicates that if the cone and sphere have the 
same slope at the boundary, the sphere is not as strong as the cone. The same 
qualitative results were also found to hold for the unrestrained cases. This 
suggests that a shallow conical cap with a dynamic buckling impulse equivalent 
to that of a spherical cap has a height larger than the sphere. 

A final comment should be made with respect to the accuracy of the present 
analysis. The two major sources of error in the analysis lie in the assumption 
that the shell behaves as a single-degree -of -freedom system and the requirement 
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that the shell deform axi symmetrically. Both of these effects become more 
important with increasing shell rise. References 1 3 and 14 present results 
for the axisymmetric dynamic snap-through buckling of a clamped spherical cap 
using five degrees of freedom. The results are shown in figure 9 as the dotted 
line. In general good correlation is obtained for the single degree of freedom 
up to A of about 7- Some disagreement between the present results and those 
of reference 1^ occurs for the cutoff point j however, this may be attributed 
to the difficulty in defining a finite cutoff point in the method used by 
reference 14. 

Results for the axisymmetric static buckling of a clamped spherical cap 
are known to be adequate up to A of about 5-5 with nonsymmetric effects pre- 
dominating for larger values of A. Since the dynamic snap-through requires a 
certain amount of time to occur, it seems reasonable to expect that axisym- 
metric behavior will hold in the dynamic case up to a higher A than for static 
snap-through. In view of these effects the results given herein are felt to be 
inadequate for A beyond about nine for conical shells and for A beyond seven 
for spherical shells. 


CONCLUDING REMARKS 

A solution has been given for the nonlinear flexural vibrations of a thin- 
walled circular cylinder subjected to harmonic excitation of a single mode. It 
was found that for certain frequencies of the driving force, only the driven 
mode responded. For other frequencies, however, the companion mode also became 
excited due to nonlinear parametric coupling between the modes. Results given 
indicate that the nonlinearity for a cylindrical shell can be either of the 
softening or hardening type depending on the aspect ratio of the vibration mode. 
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This conclusion differs from the results of Chu which indicated that the non- 
linearity was always of the hardening type. The available experimental data 
agrees reasonably well with the analysis presented here. A discussion is also 
given which relates the present study to the nonlinear vibrations of other 
axi symmetric systems. 

Results are also given for the snap-through buckling of shallow conical 
and spherical caps subjected to a uniformly distributed idealized impulse. 
Solutions are given for a variety of boundary conditions for both types of 
shells and several comparisons are made between the results. One comparison 
shows that for conical caps a simply-supported shell can resist a larger impulse 
than a clamped cone. Another comparison between spherical and conical caps 
having the sane height indicates that the sphere can resist a larger impulse 
than the cone before buckling. 
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Figure 3 .- Inf luence of large amplitudes on the vibration frequency for various values of 
Free vibrations; one mode; | * 1/2 and 2. Poisson's ratio, p * 0.3* 




Figure 4.- Influence of large amplitudes on the vibration frequency for various aspect ratios. 
Free vibrations; one mode; e = 1.0. Poisson's ratio, p = 0-3- 
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Figure 9.- Snap-through impulse versus shell parameter for axisymmetric shallow spherical 
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